Stable classical structures in dissipative quantum chaotic systems 
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We study the stability of classical structures in chaotic systems when a dissipative quantum 
evolution takes place. We consider a paradigmatic model, the quantum baker map in contact with 
a heat bath at finite temperature. We analyze the behavior of the purity, fidelity and Husimi 
distributions corresponding to initial states localized on short periodic orbits (scar functions) and 
map eigenstates. Scar functions, that have a fundamental role in the semiclassical description of 
chaotic systems, emerge as very robust against environmental perturbations. This is confirmed 
by the study of other states localized on classical structures. Also, purity and fidelity show a 
complementary behavior as decoherence measures. 
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Since the origin of quantum theory until present the 
correspondence between the classical and quantum views 
of nature has been a major source of debate. Two broad 
areas of physics have converged naturally towards its 
study. On one hand we have the semiclassical theory 
of closed systems, where the problem consists of link- 
ing the Hamiltonian with the unitary quantum evolu- 
tion as H — ► 0. For integrable systems the EBK quanti- 
zation scheme [l[ provides a precise correspondence be- 
tween quantum numbers and energy levels and between 
eigenf unctions and invariant tori. For chaotic systems 
the correspondence is not so precise but the Gutzwillcr 
trace formula Q and more recently, the short periodic or- 
bits theory Q are resources at our disposal. However, a 
complete description of chaotic eigenfunctions in terms of 
classical invariants is still lacking. On the other hand, we 
have the decoherence theory that has established a cor- 
respondence between dissipative quantum systems and 
their classical analogs These studies have enormous 
relevance nowadays not only from the theoretical, but 
also from the experimental point of view. In fact, deco- 
herence poses an unavoidable difficulty to the coherent 
manipulation of small scale systems, which in addition 
generally have a complex dynamics. Hence quantum me- 
chanics, nonlinear dynamics and decoherence are neces- 
sary ingredients in many areas of physics, involved in a 
broad range of theoretical models and implementations. 
Just a few examples are transport in cold atoms and 
BECs M, nanodevices @, microlasers @, H, @|, quantum 
dots [l(|, and chaotic scattering [ri| . 

In view of this we can ask ourselves, are there any 
classical structures embedded in the quantum realm that 
emerge as specially resistant to external perturbations? 
If so, in which way? To answer this questions, in this 
letter we investigate the stability of initial states with 
different degrees of classical information under the effects 
of chaotic evolution and decoherence. The environment 
is introduced by coupling the system to a heat bath at 



finite temperature, in a way close to actual experimental 
situations. In order to capture all the essential features 
of dissipative quantum chaotic systems without unnec- 
essary complications, we focus our study on dissipative 
quantum maps flil ] . They constitute an ideal testbed for 
semiclassical and decoherence theories. 

We find that scar functions, which are states local- 
ized along the stable and unstable manifolds of periodic 
orbits, are highly stable with respect to environmental 
perturbations. We quantify this stability by studying 
the purity and the fidelity, quantities that show com- 
plementary in order to measure decoherence. Moreover, 
by analyzing the Husimi distributions in phase space we 
conclusively show this behavior. 

We perform the evolution of the density matrix of the 
system by means of a two step operator p [H, EH ■ We 
use a composition of a unitary step given by the closed 
map B (representing the system dynamics), and a purely 
dissipative step given by the superoperator D Q , in the 
form 



S(p) = D a (B(p)) 



(1) 



In this equation B represents the unitary superoperator 
and a all the parameters of the environmental model. In 
the following we explain the construction of these two 
steps. 

We consider the baker map B on the unit torus [3| 
as the system. It is given by (q',p') = B(q,p) = 
(2q — [2q],p + [2q]/2) where (q,p) are the position and 
momentum coordinates and [a;] stands for the integer 
part of x. This transformation is an area-preserving, 
uniformly hyperbolic, piecewise- linear and invertible map 
with Lyapunov exponent A = In 2. The phase space has 
a very simple Markov partition consisting of two regions 
(q < 1/2 and q > 1/2) associated with the symbols 
and 1, for which there is a complete symbolic dynamics. 
The action of the map upon symbols can be understood 
by means of the binary expansion of the coordinates 
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(p\q) = . . . v-x ■ vqv x . . . — > = . . . v-iVq ■ vi . . . 

where q = E™o"' r( ' +1> and P = Ei=^i ^ 2 *- Then, 
a periodic orbit (PO) of period L can be represented by 
a binary string v of length L. The coordinates of the 
first trajectory point (<?o,Po) on the periodic orbit can 
be obtained explicitly in terms of the binary string as 
qo = vvv . . . — v/(2 L — l) and po = -lAiAiA . . . = 
u< /(2 L — 1), where v is the integer value of the string v 
which represents a binary number, and u* is the string 
formed by all L bits of v in reverse order. The other 
trajectory points can be easily calculated by iterations of 
the map or by cyclic shifts of v. 

When quantizing this system any state \tp) must satisfy 
periodic boundary conditions on the torus, for both the 
position and momentum representations. This amounts 
to taking (q+l\ip) = e l2jrx " (q\tp), and (p+l|^) = 
e ^x P (p\ip) , with Xqi Xp € [0, 1). This implies a Hilbcrt 
space of finite dimension N = (27rfi.) _1 . The dis- 
crete set of position and momentum eigenstates is given 
by \q 3 ) = \{]+X q )/N) (j = 0,1,.. .N - 1), and 
\Pk) = \(k + Xp )/N) (k = 0,1,. ..N - 1), labeled 
by the corresponding eigenvalues qj, pk- They are re- 
lated by a discrete Fourier transform, i.e. (pk\Qj) = 
l/VNe X p(-i(2n/N)(j + Xq )(k + Xp )) = (G%' x *). 
Throughout the paper we assume a phase space with 
anti-symmetric boundary conditions (x q — Xp — 1/2)- 
The unitary operator B that performs the closed quan- 
tum evolution is given by 
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B = G 
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N/2 



(2) 



We incorporate dissipation and thermalization to the 
quantum map by coupling it to a bath of noninteracting 
oscillators in thermal equilibrium at a temperature T. 
The degrees of freedom of the bath can be eliminated by 
means of the usual weak coupling, Markov and rotating 
wave approximations [Lsl ]. As a result we arrive at a 
Lindblad equation for the density matrix of the system p 
that can be written as a completely positive map D Q (dt) 
in the operator-sum (or Kraus) representation 



p(t + dt) = D (eiT) (dt) (p(t)) = Kp P (t)K i 



(3) 



where 



^0 = 1-5 J2 K t K » 

itf-1 

K x = y/edt (l + n(fc)) k \p k -i){ Pk \ 

k=X 
Af-1 

K 2 = ^ y/edt ra(fc) k \pk)(jPk-i\ 



(4) 



are the infinitesimal Kraus 
J2u K i K n = 1 to first order 



operators 
in dt [l7|. 



satisfying 
In these 



equations e is a system-bath coupling parameter that 
can be associated to a classical velocity dependent 
damping (at T = gives the contraction rate of the 
phase space). The population densities of the bath are 
given by n = (exp(Ai?fe/(fcsT)) — 1) , where we have 
taken E k — p\j2, and the Boltzmann constant ks = 1. 
Then, we integrate T)( £t T)(dt) numerically to obtain the 
dissipative step. 

We now briefly describe the scar function construction 
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20l | . the main tools we use to study the stability 
of classical structures. These functions are also essential 
in the semiclassical description of chaotic eigenfunctions 
3]. They are wavefunctions highly localized on the sta- 
ble and unstable manifolds of POs, and on the energy 
given by a Bohr-Sommerfeld quantization condition on 
the trajectory. We are going to use a formulation suit- 
able for maps on the torus, though they apply also to 
general flows. The first step is to define the Periodic Or- 
bit Modes (POMs) for maps, |0pqm) as a sum of coherent 
states centered at the fixed points of a given PO v, each 
one having a phase [2(| (see Figure HJa) for an exam- 
ple). Then, the scar function is obtained after applying 
a dynamical average and can be written as 



T 

ICD = E elS,A/h 008 

l=-T 



k=X 



where T stands for the number of iterations of the map up 
to the Ehrenfest time Te = IniV/A, and £„ is the classical 
action of v (see Figure 131b) for an example). It is worth 
mentioning that for T = the POMs are recovered, but 
as T is increased this function sharpens its quasi-energy 
width while extending in phase space along the stable 
and unstable manifolds of the periodic orbit. Eventually, 
it turns into a true eigenstate if the propagation time 
reaches the order of the Heisenberg time. 

In our calculations we have considered five types of 
initial states. Besides the already mentioned scar states 
and POMs, we have studied the behavior of eigenstates 
of the map having a significant overlap with scar func- 
tions. Also, we have considered states developed along 
the stable and unstable manifolds. These last two cases 
correspond to eigenstates of the momentum and posi- 
tion operators, since for the baker map the direction of 
the manifolds coincide with those of the corresponding 
axes. In addition, we have selected q and p values that 
represent points belonging to the periodic orbits shown. 
In summary, we compare states with relevant classical 
imprints (localized on POs and their manifolds) with 
strictly quantum states (map eigenstates). All results 
correspond to N = 100. 

The behavior of the coherence is studied through the 
purity, defined as P = tx(p 2 ) (tr denotes the trace). We 
find that the scar functions, the map eigenstates that 



Orbit 01 



Orbit (Kill 



.{? 0,1 -1 




Orbit 01 

I 1 I 1 I 1 I 



2 4 6 8 10 12 14 16 2 4 6 8 10 12 14 16 
Iteration (t) Iteration (t) 



Orbit 00 11 
1 I 1 I 1 I 1 I 




6 8 10 12 14 16 
Iteration (t) 



Figure 1: (Color online) Purity P (logarithmic scale) as a 
function of time (measured in units of map steps) . In the left 
column we show results for the orbit v — 01, while in the right 
one for v = 0011. Upper panels correspond to e — 0.001, 
and lower ones to e — 0.01. Orange (lighter) lines corre- 
spond to T = 1 and blue (darker) lines to T = 1000. Solid 
lines stand for scar functions, long-dashed lines for eigenstates 
of the map, short-dashed for POMs, dot-dashed for position 
eigenstates and dotted for momentum eigenstates. N = 100. 



have a maximum overlap with them, and the POMs, all 
behave in a similar way. However, the purity decay is 
faster in the case of the eigenstates when compared to 
that of the corresponding scar functions. On the other 
hand the position eigenstates are the less and the mo- 
mentum eigenstates the most decoherent ones. We show 
results for two representative orbits in the left and right 
columns of Figure [1] (they are denoted by v = 01 and 
v = 0011, respectively). In obtaining these results we 
have considered two different couplings with the envi- 
ronment (e = 0.001 and e = 0.01), and two different 
values of the temperature (T = 1 and T = 1000). They 
convalidate the hypothesis that entropy production, then 
entanglement with the environment, is governed by the 
classical unstable behavior of the system [4(. The more 
localized on classically unstable regions the initial distri- 
bution is, the fastest its entropy production results. 

In order to have a complementary stability criterion we 
have explored the fidelity calculated as F = yj (ip\p\ip). 
This quantity measures the degree of overlap between 
the pure initial state \ip) and the evolved density matrix 
p. This time we find that although the scar functions 
and the corresponding map eigenstates have a similar be- 
havior, the latter are more stable than the former when 
the temperature is low. For higher values the situation 
changes and the scar functions become more stable than 
the eigenf unctions. The coupling with the environment 
seems to play a similar role than temperature, since the 
decay rate difference for low values of T shrinks as e 
grows. Finally, POMs fidelity decays faster, and the po- 



Figure 2: Fidelity F (logarithmic scale) as a function of time 
(measured in units of map steps). We show results for the 
same orbits, use the same parameter values, and colors and 
patterns criterion as in Fig. [1] 



sition and momentum eigenstates have the greatest decay 
rate but with big oscillations, whose periodicity is given 
by the dynamics. Hence, for low temperatures the stabil- 
ity is governed by the quantum mechanics of the system. 
Initial states that are closer (greater overlap) to the map 
eigenf unctions turn out to be the most stable ones. In 
this sense, scar functions are more stable than POMs, 
which in turn are more stable than position and momen- 
tum eigenstates. When the effects of the environment 
become stronger scar functions emerge again as remark- 
ably stable, having the lowest decay rate of all the initial 
states we have considered. This is shown in Fig. [21 where 
we have used the same orbits, parameter values, colors 
and patterns as those chosen for Fig. [TJ 

The picture is completed by means of the Husimi dis- 
tributions obtained at different times. They can be seen 
in Figure [31 where the evolution of the POMs, the scar 
functions and the map eigenstates for the orbit v = 0011 
can be found. The high stability of the scar functions 
is undoubtedly demonstrated. When already the POMs 
shape cannot be distinguished anymore they still display 
a high level of detail, showing their characteristic local- 
ization along the manifolds of the corresponding orbit. 
Moreover, the eigenstate rapidly loses its original details 
and seems to converge to the scar function. 

In conclusion, we have verified a great degree of com- 
plementarity between purity and fidelity as decoherence 
measures. The loss of purity, or entropy production, is 
governed by the classical instability of the system, while 
the fidelity decay for low temperatures (and weak cou- 
pling with the bath) is dominated by the overlap of the 
initial state with the map eigenfunctions. When tem- 
perature rises the situation changes and the localiza- 
tion on classical structures plays again a fundamental 
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Figure 3: Husimi distributions for the POM and scar function 
associated to orbit v = 0011 ((a) upper and (b) middle panels, 
respectively), and corresponding map eigenstate ((c) lower 
panel). We have used e = 0.01 and T = 1000. Time steps go 
from to 7 from left to right and top to bottom. The grayscale 
goes from white (minimum probability) to black (maximum 
probability) . 

role. By analyzing these quantities we have found clas- 
sical structures in quantum mechanics that are specially 
stable when subjected to environmental perturbations. 
They are the scar functions, which are associated to peri- 
odic orbits and the stable and unstable manifolds in their 
vicinity. Husimi distributions in phase space allowed us 
to clarify this picture. We have conclusively shown that 
scar functions keep their shape virtually intact at times 
when the POMs have already lost all their characteristic 
features. Hence, more than just localization on the pe- 
riodic points is needed to provide stability. Finally, by 
exploiting the simplicity of the manifolds of our map we 
could verify that any kind of localization on them is not 
enough to guarantee robustness against external pertur- 



bations, but the one provided by scar functions. 

Remarkably, the purity and fidelity loss of the map 
eigenstates (in this last case with the exception of low 
temperatures) is generally faster than that of the scar 
functions. Then, these latter represent the stable classi- 
cal skeleton of the map eigenstates against environmental 
perturbations. It is worth underlining that these classi- 
cal structures which survive for longer times are also the 
main classical ingredient needed to construct the eigen- 
states, according to the short periodic orbit theory. 

Support by CONICET is gratefully acknowledged. 
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